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ABSTRACT 

The nearby dwarf starburst galaxy NGC 5253 hosts a deeply embedded radio- 
infrared supernebula excited by thousands of O stars. We have observed this source 
in the 10.5/zm line of 5 +3 at 3.8 kms" 1 spectral and 1.4" spatial resolution, using the 
high resolution spectrometer TEXES on the IRTF. The line profile cannot be fit well 
by a single Gaussian. The best simple fit describes the gas with two Gaussians, one 
near the galactic velocity with FWHM 33.6 kms -1 and another of similiar strength and 
FWHM 94 kms -1 centered ~ 20 kms -1 to the blue. This suggests a model for the 
supernebula in which gas flows towards us out of the molecular cloud, as in a "blister" 



or "champagne flow" or in the Hll regions modelled by Zhu (2006). 



Subject headings: galaxies: star clusters-galaxies: starburst 



1. Introduction 

Super star clusters are a common mode of star formation in starburst galaxies. These star 
clusters form deeply embedded in molecular clouds, and the young stars disperse the remaining 
cloud material to become visible clusters, while contributing to the ionization, mechanical energy, 
and metal enrichment of the starburst. Through the course of their short lifetimes, O and B stars 
have the potential to massively change their environments through the action of stellar winds and 
supernovae, which could have the effect of either halting further star formation, or inducing a new 
round of cluster formation. What are the observed effects of a young super star cluster on its 
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immediate environment? To address this question, we are carrying out observations of the gas 
kinematics of the supernebulae these star clusters excite. This paper examines the gas kinematics 
in one of the closest such supernebulae, which is also the most luminous known. 

The dwarf galaxy NGC 5253 has created towards its center dozens of young, bright super star 
clusters (SSCs) with ages of a few to tens of Myr ( Meurer et al |1995 Calzetti et al.||1997 Harris et 



aL"|200 4 ) . In the same area as the young clusters is an infalling streamer of CO ( Turner et al.f l997 



Meier et al.|2002[ )) and HI ( |Kobulnicky Skillman|2008| |L6pez-Sanchez et al.|2007[ ) and coincident 



with it, a filament of emission first detected in [OIII] by Graham (1981) and lately interpreted by 
Zastrow et al. (2011) as an "ionization cone". It is possible that the infalling molecular stream 



encouraged formation of the clusters, and that the cone is the outside of the same cold gas streamer, 
lit up by the ionization. This illustrates how clusters can influence the galaxy well beyond their 
immediate vicinity. 



Infrared observations ( Turner et al.|2003 Alonso-Herrero et al.|2004 ) of the core of the starburst 
region find a deeply embedded star cluster located 0.3 — 0.4" (5-6 pc at 3.8 Mpc) to the northwest 
of the brightest optical cluster (cluster "5" of Calzetti et al. ( |1997 ).) The embedded cluster is not 
detected at wavelengths shorter than 1.9 ^m due to extinction, which is estimated to be Ay ~ 17, 



and appears to originate within the source (Turner et al. 2003 Martin-Hernandez et al. 2005) 



The embedded cluster is coincident with bright Brackett line emission associated with a strong 
radio source, the "supernebula" ( |Turner et aLlp98"l [2000} |Turner fc BecF|[2004l ) . The supernebula 
requires the ionization equivalent of 1200 07 stars to excite the dense, parsec-sized core of the 



HII region (Turner & Beck 2004; Mohan et al. 2001). In addition to the dense core, there is a 



more extended halo to the H n region, requiring the excitation of 2000 07 stars within the central 



5 pc, and up to 7000 O stars within a radius of 20 pc (Turner & Beck 2004). The extinction, 
the high inferred density of several thousand per cm 3 for the Hn region (Turner & Beck 



2004 



Monreal-Ibero et al. 2010) and the young spectroscopic age for the cluster (Alonso-Herrero et al. 



2004) suggest that this is one of the youngest known superclusters, still in its embedded phase and 



possibly still in the process of formation. 

The current star forming episode in NGC 5253 is dominated by the "supernebula", which is 
only ~1 pc in radius. A luminosity in young stars of ~ 10 9 L and ionization rate of ~ 4-7 x 10 52 s _1 
within a 2-10 pc region require the presence of a stellar cluster containing thousands of O stars and 
a dense, dusty, compact H n region. This is not a static configuration: hot ionized gas will expand 
under its high pressure, and young stars, especially O stars, are active and drive powerful winds 
and outflows. What effect is this young, forming super star cluster having on its surroundings? 



Are stellar winds evident? If there are O stars present, especially evolved O stars (Alonso-Herrero 



et al. 2004), then why is this HII region still embedded? What can the gas kinematics be in such 
an extreme Hn region? 



The kinematics of the supernebula have been observed in the near-infrared Brackett recom- 
bination lines of HI (Turner et al. 2003) . The FWHM of the near-infrared lines, measured with 
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12 km s 1 resolution, was 76 ± 2.5 km s 1 over the central few arcseconds surrounding the superneb- 
ula. These velocities are not much higher than the FWHM of some individual Galactic compact 



H II regions excited by one or a few O stars ( Beck 2008 ) . But compared to a normal H n region 



the supernebula has thousands of times more stellar mass in roughly the same volume. This led 



Turner et al. (2003) to suggest that the expansion velocities may actually be less than the escape 
velocity for much of the cluster, so the supernebula could be gravitationally bound; a giant photo- 
sphere, rather than a conventional H n region. This would correspond to the "Type III" cluster of 



Kroupa & Boily (2002), in which the escape velocity exceeds the sound speed. On somewhat larger 
(60-100 pc) scales, gas motions have also been observed through the H53a radio recombination line 
( |Rodriguez-Rico et al |2007[ ) , and in red lines of Ha, [N II] and [S II] ( |Monreal-Ibero et al.|2(5To| ) . The 
H53a line was measured with 44kms~ 1 resolution and had an observed FWHM ~ 100 ± 12kms _1 



over the inner 4" (~60 pc, see Figure 5 of Rodriguez-Rico et al (2007)). Rodriguez-Rico et al 
(2007) found a gradient of 10 km s _1 arcsce -1 over this region, which they interpreted as rotation. 



Monreal-Ibero et al. (2010) also find a gradient, over scales similar to the H53a image, but it is a 
gradient in differential velocity between the forbidden lines [Nil] and [S II] and Ha. Neither the 
H53a nor the red lines are ideal for probing spatial and kinematic effects in the the supernebula; 
for the RRL the beam is too big (2") to distinguish motions within the nebula while extinction 
within the supernebula prevents observation in the optical. 

The strong mid-infrared emission lines of the metal ions are useful probes of ionized gas in 
obscured regions. They are relatively free of extinction compared to visual lines and the low thermal 
broadening due to the ionic weight means that they are much better at probing kinematics than 



Zhu et al. 


2008; 


Jaffe et al. 


2003) 



of 10 4 K the thermal width of a hydrogen line will be 21kms _1 FWHM, but the thermal width 
V m of line emitted by a metal of mass m m will be smaller than that of a hydrogen line by the 
ratio V m /Vh = (mh/mm) 1 ^ 2 . For NGC 5253, which is a high-excitation source, the strongest line 



at mid-infrared wavelengths that can be observed from the ground (Beck et al. 1996) is [SIV] 
10.5 /um. This promises a factor of 5.66 reduction in the thermal component of the line profile. 
We accordingly observed [SIV] 10.5 /im from the supernebula in NGC 5253 with high spectral and 
~ 1" spatial resolution, with the aim of separating the thermal component of the line width and 
determining the true bulk motions of the gas. 



2. Observations; the TEXES [S IV] Data Cube 

NGC 5253 was observed on the night of 4 June 2010 with TEXES, the Texas Echelon Cross 
Echelle Spectrograph ( |Lacy et al.]|2002| ) spectrometer, on the NASA IRTF on Mauna Kea. TEXES 
is a sensitive spectrometer for the 5-25 ^m region, with three resolution modes: these data were 
obtained in the high resolution mode, which gives spectral resolving power R ~ 80, 000. The slit 
was 29 pixels long and the plate scale 0.36" x 0.95 km s _1 per pixel. Diffraction and seeing limit 
the final spatial resolution to 1.4" and the spectral resolution to ~ 3.8 km s -1 . We offset from the 
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guide star to the position of the radio supernebula (Turner &; Beck 2004). The slit, which was 
oriented NS, was stepped across the galaxy in 1" increments, first west, then back to the center to 
check the pointing, which was excellent, and then east. At each position the galaxy was nodded 
12" north, taking the emission region off the 9" slit. The beams were subtracted and the spectra 
were combined to create a data cube. 

Wavelength calibration is done by reference to Earth atmospheric lines seen in the spectrum 
of the source and checked against the spectrum of the asteroid Vesta, which was also used as the 
divisor to remove the atmospheric transmission. There are no strong lines close to the wavelength 
of [SIV] in NGC 5253, which from the redshift of ~ 400 kms -1 and the [SIV] rest wavelength of 
951.43 cm -1 is predicted to be 950.16 cm -1 . Velocities quoted are heliocentric. 

The total flux in the [S IV] map (the zeroth moment of the data cube) is 6 x 10 -12 erg s _1 cm -2 , 



50% higher than the previous result oflBeck et al. ( 1996 ). This may reflect a more thorough mapping 



of the [S IV] emitting region or an offset in the calibration. This flux agrees, within the calibration 
uncertainties, with the total [S IV] emission seen by the Spitzer IRS in a much larger slit. 



3. Spatial Distribution of the Gas 

The TEXES data give us effectively a spectrum at each point in the sky. The data were 
obtained with the slit oriented NS, so by collapsing the spectra along the slit we obtain 1-d spectra 
showing the line profile in each EW position. The line was detected in each position over the total 
spatial extent of 5", but the strongest position is more than 5 times stronger than any other. The 
spectra are displayed in Figure 1. The strongest position is 1" east of the point arrived at by offset 
from the guide star to the radio source. In light of the spatial resolution, beam size, and usual 
offset accuracy this is not significant and we identify the peak of the [S IV] emission with the radio 



supernebula at : R.A.= 13 h 39 m 55 s .96, decl. = -31°38'24.38" ( |Turner fc Beck ||2004| ). In Figure 
2 we show the [S IV] data cube in the form of a grid of spectra, obtained by binning the data into 
1" x 1" cells. The grid steps are comparable to the diffraction limited beam size of the [S IV] data. 

The [SIV] is concentrated but it is not point-like. In Figure 3 we show the spatial distribution 
of the [SIV] emission integrated over the whole velocity range (the zeroth moment of the data 
cube). The [SIV] emission is symmetric and only slightly extended EW, but in declination it is 
clearly extended and asymmetrically so: the emission falls off sharply north of the peak but extends 
3" towards the south. This is seen along individual NS slits and is not a result of the observing 
procedure or the data reduction. The extension of the [S IV] line emission south of the supernebula 



agrees with the high spatial resolution radio continuum map of Turner Sz Beck (2004). Their 
image shows a "filamentary arc" of free-free radio emission extending south of the main source 
and appearing to curve around the second nuclear star cluster (the less obscured cluster which is 
optically visible but not a strong radio source); they suggest that it may be a blister feature or 
cometary Hll region. 
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4. The [S IV] Line Profile 

It is clear from Figure 1 that the [SIV] line profile is not simple and symmetric. There is 
obvious red-blue asymmetry around the peak, with the blue side extending to higher velocity than 
the red, and the whole line profile appears distinctly non-Gaussian. Given the complex line shape, 
and that there are no models for the gas kinematics of a supernebula or for what the resulting line 
profile 'should' look like, caution is necessary in trying to fit the line profile. While it is possible 
to obtain a fit of great formal goodness by using many components, it may not be physically 
meaningful. Here we concentrate on the simplest fits. They are not formally perfect, but can still 
reveal a great deal about the supernebula. 

Since we do not a priori know the kinematics of gas in the supernebula, the first step is 
to test what aspects can and cannot be described by the simplest possible fit, which is a single 
Gaussian. We used the IDL procedure LMFIT to search for the fits. The best single Gaussian 
fit to the [SIV] line, and the formal uncertainties to the fit, are: center at +377 ± 0.5 km s" 1 and 
a = 27.8 ± 0.6 km s -1 , or FWHM of 65kms _1 . This fit and the residuals are shown in Figure 4. 
The residuals in Figure 5 clearly show that two features of the [S IV] profile are not well fit by a 
single Gaussian: a narrow component about lOkms" 1 red of the main peak, and excess emission 
on the blue side extending out to -100 kms -1 . 

Two Gaussians, shown in Figure 5, produce a much better formal fit. One Gaussian is centered 
very close to the galactic velocity at 391 ± 0.8 km s -1 and has a = 14.1 ± 0.7 km s -1 or FHWM 
33.6 km s" 1 . The second, of almost equal peak strength, is offset 17±6kms" 1 to the blue and has 
a = 39.25 ± 1.7kms _1 or FWHM 94kms~ 1 . The formal reduced x 2 01 the single Gaussian fit is 
2.04 and of the two Gaussian fit, 1.29. The absolute values of x 2 depend on the noise level and may 
have systematic errors; the relative values of x 2 f° r the fits are more significant, and agree with the 
clearly superior appearance of the two Gaussian fit and its residuals. We use these fits as guides 
for our discussion of the nebular structure; as more information on the structure and kinematics of 
the supernebula, the cluster and its environment is recorded, a more realistic model can emerge. 

Gaussian line shapes are produced by random motions. The gas in the supernebula is certainly 
undergoing at least two kinds of random motions: thermal, which depend on T e , and virial turbu- 
lence which depends on the gravitational potential of the star cluster mass. There may be other 
sources of random motions and line broadening present. For example, the intense star formation 
that created the embedded cluster was presumably accompanied by stellar winds and outflows that 
could have generated significant turbulence. But these other sources, if they exist, are not yet 
quantifiable. We therefore discuss here only the gravitational and thermal effects. 
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4.1. Gravitational Effects 

How do the velocities from the fit compare to estimates of the gravitational energy of the 
system? From the and total luminosity of the supernebula and a Kroupa IMF extending to 
O.IMq, the embedded cluster has total mass ~ 3 x 1O 5 M . Using the virial velocity dispersion 
ay = (GMci/311) 1 / 2 and assuming that the Gaussian component of the line emission originates in 
the same compact region as the free-free continuum, which we take from the image of Turner & 



Beck (2004) to be R ~ lpc, gives a r ~20kms -1 or FWHM for the bulk virial motions ~ 47 kms -1 . 
This rough calculation shows that the bulk of the gas has velocities consistent with other estimates 
of the total cluster mass; there is no need for a truncated or top-heavy IMF to explain the M/L 



ratio. 



4.2. Thermal Broadening and the Profiles of [S IV] and Brackett Lines 

Any nebular line will have a thermal width which depends on the ionic weight and T e . A 
motivation to observe the [S IV] line is that sulfur's weight makes it less susceptible to thermal 
broadening than are hydrogen lines, so it is possible to see the gas motions in more detail. How does 



the [SIV] line compare to the Brackett line profiles of Turner et al. (2003)? We will concentrate 
on the Bra line because of the longer wavelength and lesser susceptibility to extinction, since 
A2^ m ~ 2 mag, and we sum the Bra over 1" to match the [SIV] slit. The best single Gaussian 
fit to the Bra line is centered at 377 km s -1 and has a = 31.5 kms , FWHM = 74 kms -1 . The 
best two Gaussian fit has a narrow component, a = 22.7 kms -1 , centered at 391 kms -1 and a blue 
component with a = 49.1 kms -1 centered at 368 kms -1 . The best single Gaussian fit and residuals 
are shown in Figure 6 and the best two Gaussian fit and residuals in Figure 7. For Br a the single 
Gaussian and two Gaussian fits have the same reduced x 2 ~ 0.96. This is a striking contrast to the 
[S IV] line, which could clearly not be fit well by one Gaussian. 

The [SIV] and Bra lines are expected to probe the same body of gas; while S +3 and H + have 
different ionization potentials they are expected to co-exist in a nebula ionized by hot stars because 
almost all the sulfur will be in S +3 . This is supported by the ionic abundance of S +3 derived from 
the [S IV] line strength, which is at least as high as the elemental abundance of sulfur in NGC 5253 



(assumed to be ~0.29 solar, Kobulnicky et al. (1997)'s value for the galactic metallicity). So the 



lines should resemble each other except in so far as the Brackett lines have lower spectral resolution, 
greater thermal broadening and more sensitivity to extinction. That two Gaussian fits match both 
the [SIV] and Bra profiles show a basic agreement of the two lines. There is more evidence for 
the blue wing in [SIV] than in Bra, even degrading the [SIV] resolution to 32 kms -1 (the sum 
in quadrature of the instrumental resolution and the thermal width of H at 2 x 10 K, see next 
section), leaves more blue emission than is seen in the Bra. But the Bra line is affected by a 
negative feature 70 kms -1 blue of the peak, perhaps imperfectly cancelled atmospheric absorption, 
which could disguise a blue excess. We would need better Br a data to be sure that the line profiles 
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are consistent. 



4-2.1. T e in the Supernebula 

One can derive the thermal broadening of a spectral line, and thus T e of a nebula, by comparing 
the line widths of ions of different weights. We use for this the single Gaussian fits; although the 
single Gaussian is not an optimum fit for the [S IV] it is most appropriate for the current Br a 
line data. If higher resolution Brackett spectra become available, it may be possible to fit the line 
better and derive T e more accurately. 

If we compare the single Gaussian line fits for hydrogen and sulfur and assume that the widths 
are due to the convolution of a thermal width, which is 5.66 times greater for hydrogen than for 
sulfur, with a turbulent or other random width which is the same for both species, the resulting 
simultaneous equations can be solved to give a thermal width of 5.8 km s -1 for sulfur and 32.5 km s _1 
for hydrogen. This is significantly higher than the 21 kms^ 1 FWHM velocity broadening associated 
with a 'normal' H n region temperature of 10 4 K and requires T e of 2.5 x 10 4 K. This T e value is quite 
realistic for the supernebula; T e depends on ionization parameter U, and the small size and intense 
radiation of the supernebula create a very high ionization parameter, log U ~ —0.1. Donahue & 



Shull (1991) calculate cooling in ionized clouds and predict a minimum nebular temperature of 



2 x 10 4 K for this U and the metallicity of NGC 5253. 



4.3. Modelling the Blue Wing 

Both the single Gaussian and two Gaussian fits indicate excess blue emission. Blue wings 
are often observed in the spectra of embedded, highly obscured stars and clusters. An obvious 
model for this profile in a Hn region is that the blue wings are due to supression of the red side 
by extinction from absorbing dust mixed with the emitting gas in a symmetric expanding flow. 
Arsrv] ~ 0.7Ak, which would give Arj/y] of 1.4 magnitudes: a significant amount of obscuration, 
worse than the Bra suffers. But it should be noted that the silicate feature is weak in NGC 5253 
and this figure is an upper limit, possibly quite a high one, for the [SIV] obscuration. But while 
there is ample extinction in the supernebula ( Turner et al.|2003 found A^os^m to be 0.8 mag) the 



model is not very physically plausible. First is a lifetime problem: if the supernebula were really 
expanding at ~ 100 km s -1 , as the linewidths would imply in this model, it would have achieved the 
radio size of 0.7 pc in ~ 6 x 10 3 years and the 1" of the [SIV] resolution in only 1 x 10 5 years old. 
Yet observations of the star cluster and the evidence for possible WR activity in the area suggest 
that it is ~ 3Myr old ( Alonso-Herrero et al.|[2004 ). The second argument against this model is the 



line shape itself. If the absorbing material in the expanding cloud is uniformly distributed, then 
extinction sufficient to reduce severely the red-shifted emission would be expected to distort and 
shift the line center relative to other lines with different obscuration. But the peak velocity of the 
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[SIV] agrees with the Brackett lines (Turner et al. 2003), which agree with each other, in spite of 
their different extinctions. 

Another possible explanation of the blue wing is that the emitting region contains a secondary 
source or sub-cluster blue-shifted relative to the main cluster. If so, the line shape should vary with 
position, except in the unlikely case that the secondary source is exactly in the same line of sight 
as the main cluster. We analyzed the spatial dependance of the line shapes by interpolating onto a 
grid 0.36 x 0.33" pixels (close to the 0.36" size of a detector pixel). In Figure 8 we show the spectra 
obtained in the peak pixel and in a 3 x 3 pixel box (minus the peak pixel), a 5 x 5 minus the inner 
3x3 pixels, 7x7 minus the inner 5x5 and 9x9 minus the inner 7x7 pixel boxes surrounding 
it. To the limits of the signal-to-noise, no variation is apparent. We believe that while a secondary 
blue-shifted source cannot be ruled out, there is at present no evidence for it. It should be noted 



that high resolution radio maps (Turner et al. 2000) do find another small source 0.3" from the 



main cluster, but it is a much smaller fraction of the total radio emission than the blue wing is of 
the [SIV] flux. 



.3.1. Instrinsically Asymmetric Motions-A Blue Flow'; 



The supernebula line profile is not unique: Henry et al.| (2007) fit Brackett lines in the super star 
cluster nebulae in He 2-10 with narrow (60- 71km s" 1 ) plus broad (200 - 300 km s" 1 ) profiles; the 
broad profiles are blueshifted with respect to the narrow components, which they attribute to flows 
in an aging starburst. Gilbert & Graham (2007) report Br7 lines with FWHM 60 — llOkms -1 and 



non-Gaussian wings up to ~ 200 km s -1 wide for 17 embedded cluster sources in the Antennae 
galaxies, also in regions that are ~ 100 pc extent, several times larger than those we consider here. 
Since supernebulae are in many respects (e.g. radio spectrum, evolutionary stage) scaled-up version 



of Galactic Ultra-Compact Hil regions, and Beck (2008) showed that the HI recombination line 



profiles of embedded clusters resemble the profiles obtained by superposing the HI emission of 
many ultra-compactH n regions, it is natural to look at these H n regions-which can be studied in 
great spatial detail-for clues to the possible gas motions. 

The infrared emission lines of UCHll regions frequently have a narrow line center and wide 
low level flux, over velocity ranges like those of the [SIV] lines, and are asymmetric, in the same 



sense (blue excess) as we see in NGC 5253. Persson et al (1984) present Bra spectra of embedded 
young stellar objects which have deconvolved FWHM 50 — 235kms _1 and, in some cases, broad 



blue wings. Kurtz (2005) shows radio recombination lines of similar FWHM from hyper-compact 



and ultra-compact Hn regions. High resolution spectral mapping (Zhu 2006) and modeling has 



shown that in these cases the ionizing star(s) lie on the near side of the natal molecular cloud and 
that gas flowing along the cloud surface and away from the cloud creates the blue wing. 



We think this is the most likely picture for the supernebula and the most natural explanation 
of the line shapes. It explains why the core of the line is almost symmetric; that gas has not 
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entered the flow zone and is behaving like a normal H n region. We note that although Wolf-Rayet 
winds are not thought to turn on until ~ 3 Myr, Zhu et al. (2008) found that most observed 



ultra-compact H n regions show evidence of the influence of stellar winds. 



5. Discussion and Conclusions 



The supernebula in NGC 5253 is the most compact embedded super star cluster known, one 
of the most obscured, and one of the very brightest. It is also probably the youngest, as shown by 
its small size, high gas density (about 3 — 4 x 10 4 cm -3 , Turner & Beck 2004), and the very high 



ionization state of the gas (Crowther et al (1999) arrived at an age of 2-3 Myr from analysis of 



the radiation field, agreeing with 3 Myr Alonso-Herrero et al. (2004) found from the near-infrared 



continuum). It is in a complex region of many star clusters and molecular gas. As it ages, emerges 
from its embedding material, and produces supernovae it will affect its surroundings, and indeed 
the whole galaxy, profoundly. The kinematic information gained from the [SIV] line has let us 
probe the gas motions in more detail than ever and can point out how the supernebula may evolve. 

The [S IV] spectrum we present here measures the ionized gas excited by an embedded star 
cluster with the highest spectral resolution yet achieved. The [S IV] line shapes are quite constant 
throughout the [S IV] emission region, with no sign of the apparent gradient seen on larger scales 



by the radio recombination lines (Rodriguez-Rico et al 2007), nor those seen in the red forbidden 



lines of [S II] and [Nil] ( |Monreal-Ibero et al.|2010 ). The line profile is not a simple Gaussian, but is 
asymmetric about the peak with a blue wing. We look at the simplest possible fits as starting points 
for discussion of the kinematics. The single Gaussian which is closest to fitting the line has FWHM 
65kms _1 ; this fit leaves excess blue emission and an additional narrow (22kms _1 ) redshifted peak. 
A much better formal fit is obtained with two Gaussians, one of FWHM 33.6kms _1 at galactic 
velocity and another with FWHM 94kms _1 offset 17 ± 6kms _1 to the blue. There are no clear 
features in the residuals of this fit. 

The best formal fit leads to a model in which the gas may be viewed as either still confined to 
the cluster (the main peak), or having entered the flow away from the cluster (the blue component). 
In reality these are not two distinct bodies of gas, but continuous. In this model the supernebula's 
kinematics are like those of a blister H n region, where the ionized gas flows out of the cloud. Note 



that when Turner & Beck (2004) spatially resolved the ionized gas of the supernebulae, they found 



an arc or cometary shape, which again resembles a blister H n region. As more information on the 
supernebula and on the surrounding clouds becomes known, it should be possible to refine the fits 
according to what is physically most realistic, as well as formally good. 

This blister or pressure-driven flow picture is an attractive one for supernebula. Because the 
ionized gas does not simply expand, there is no life-time problem. The overpressure of the ionized 
gas is relieved by the outflow, while the supernebula maintains its small size, and this can continue 
as long as the gas flowing out is replenished by ionization of molecular gas from the cloud. Other 
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extragalactic supernebula should be observed via the infrared metal lines to determine if they show- 
similar line profiles at high spectral resolution. 

TEXES observations at the IRTF were supported by NSF AST-0607312 and by AST-0708074 
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Fig. 1. — Top: The spectrum in the peak position averaged over 3" along the slit. Intensity units 
are erg(s cm -1 cm 2 sr) _1 ; x axis is heliocentric velocity. Bottom: Spectra in the weaker positions 
averaged over 3" along the slit. The spectra are, from top to bottom, 1" east, 1", 2", and 3" west 
of the peak. Spectra are offset vertically by steps of 0.015. Units as in previous figure. 



-14- 



.1 - 

05 - 

06 ~ T 
04 - 
02 ~ T 

-i i i i 1 ICi i i 1 i i i i 1 i 1 T i 1 i i i 










05 ~ r 

06 - 
04 - 

; , , , i 


St 








oa - 

06 ~ r 

04 r 

02 F- 

; , I , , ,1 


:, , , , 1 , . , , 1 , , , , 1 . , , , 1 , , , 








08 '- 
06 - 
04 - 

02 - 

; , , i , ! , ! 











100 200 300 400 500 600 200 300 400 500 600 200 300 400 500 600 200 300 400 500 600 200 300 400 500 600 



Fig. 2. — The spectrum of [SIV] at each point on a 1" x 1.1" grid. North is up and East left. The 
data were binned by 3 pixels along the slit and 2 pixels in velocity. Units as in Figure 1. 



Fig. 3. — [SIV] line intensity integrated over +275 — 1-465 km s 1 , showing the spatial distribution 
of the line. North is up, East is left. Contours are integer multiples of 1.6 x 10~ 2 erg(s cm 2 sr) _1 . 
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Fig. 4. — The [SIV] line observed in the strongest position, the best fitting single Gaussian, and 
the residuals. 
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Fig. 5. — The [SIV] line observed in the strongest position, the best fitting sum of two Gaussians, 
and the residuals. 
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Fig. 6. — The Bra line observed with 12kms _1 resolution. The data is from Turner et al. (2003) 
but taken over a 1.2" slit. Also shown are the best single Gaussian fit to the line and the residuals 
of the fit. Y-axis units are counts and velocity is LSR. 
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Fig. 7. — The Bra line of Turner et al. (2003), the best two Gaussian fit and the residuals to the fit. 
For this line, in contrast to the [S IV], the two Gaussian fit is not clearly better than one Gaussian. 
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Fig. 8. — The [SIV] line in boxes or annuli around the peak. The boxes are 3x3, 5x5, 7x7 and 
9x9 pixels. The smaller boxes are subtracted from the larger ones, i.e. the 5x5 pixel spectrum 
does not contain contributions from the inner 3x3 pixel box. Units are as in Figure 1. 



